Local control of ultrafast dynamics of magnetic nanoparticles 



oo 

o 

O 

(N 

O 

Q 

(N 



5-H 



-a 
c 

o 
o 



> 

00 

in 
o 

(N 
00 

o 



X 



A. Sukhov 1,2 and J. Berakdar 2 

1 Max-Planck-lnstitut fiir Mikrostmkturphysik, Weinberg 2, D-06120 Halle/Saale, Germany 
2 Institut fur Physik, Martin-Luther-Universitat Halle-Wittenberg, Heinrich-Damerow-Str. 4, 06120 Halle, Germany 

(Dated: December 2, 2008) 

Using the local control theory we derive analytical expressions for magnetic field pulses that steer the mag- 
netization of a monodomain magnetic nanoparticle to a predefined state. Finite-temperature full numerical 
simulations confirm the analytical results and show that a magnetization switching or freezing is achievable 
within few precessional periods and that the scheme is exploitable for fast thermal switching. 

PACS numbers: 75.10.Hk, 82.50.Nd, 75.60Jk , 75.40.Mg, 87.19.lr 



Introduction.- A fast magnetization reversal of magnetic 
nanoparticles is of a key importance for the realization of 
high-rate magnetic recording JjJ, |2y. Several techniques are 
currently envisaged for the magnetization switching such as 
the laser-induced spin dynamics [3] based on the inverse Fara- 
day effect 

Si, 

the reversal triggered by external static or 
alternating magnetic fields SHSSllSliiEl or by a spin- 
torque acting on the magnetization due to a passing spin polar- 
: lil3lll4ll . Transverse magnetic field pulses 



16, 17 



ism, 



ized electric current 
are also efficient for a swift reversal 

and if finely tuned in duration J2, 12. lh can even lead to a 
quasi-ballistic switching. A further fundamental issue, ad- 
dressed here is how to steer the magnetic dynamics to a de- 
sirable state by external fields. Generally, a number of control 
schemes have been established mainly in quantum chemistry 



JH m m 

theory (LCT) 



Particularly interesting is the local control 



24 I25H in which the control fields are con- 



structed from the response of the system offering thus a physi- 
cal interpretation of the control mechanism. We adopt the idea 
of LCT to steer the magnetization dynamics of nanoparticle by 
transverse magnetic pulses. We obtain transparent analytical 
expressions for the control pulses that allow a fast switching or 
a quasi "freezing" at a predefined magnetization state. For the 
scheme to be applicable, the field durations have to be shorter 
than the field-free precessional period but no special pulse- 
duration tuning is required; the field strengths are to be deter- 
mined according to the analytical expressions provided here. 
In our control strategy the magnetization dynamics proceeds 
via sudden impulsive kicks guiding the magnetization towards 
a predefined direction; the pulses are intervened by field-free 
magnetization precessions and relaxation. A similar mech- 



anism has recently been realized experimentally [12] using 
spin-polarized picosecond current pulses resulting in a spin- 
transfer-torque-driven stroboscopic dynamics. The robustness 
of the predictions we demonstrate with finite-temperature full 
numerical calculations and for different types of anisotropy 
fields. We confirm the analytical results and uncover the po- 
tential of this scheme for fast thermal switching that can be 
the basis for fast thermal sensors. 

Theory.- We consider a nanoparticle with a size such that it 
displays a long-range magnetic order and is in a single domain 
remanent state. Examples are FesoPtso J3, 13l or Fe7oPt3o 
J2, 01 nanoparticles which possess respectively a uniaxial 



or a cubic anisotropy. Following the Landau-Lifshitz-Gilbert 
(LLG) approach we model the dynamics of the magnetization 
direction by the classical evolution of a unit vector S. The 
particle's magnetic moment at saturation /ig is assumed time- 
invariant. The system energy derives from H = Ha + Kf, 
where Ha and Hf = — S ■ bo(t) stand respectively for the 
anisotropy and the Zeeman energy of S in the external field 
ho(t). For a particular type of anisotropy described by /a(S) 
we write Ha = —DJa{S) with D being the anisotropy 
constant. S(t) develops according to LLG equation B28I1 

where 



as 



as 

01 



B e (f) + a(S x B e (t)) 



(l+a 2 )^ 

B e (t) = -l/(ps)9H/dS is the effective field, 7 is the gy- 
romagnetic ratio and a is the Gilbert damping parameter. In 
spherical coordinates where the z axis along is the easy axis 
we specify S by the azimuthal (</>) and polar (9) angles and 
cast the LLG equation as lul 12911 
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Hereafter the time is measured in units of the field-free preces- 
sional period T pl0C and the energy H in units of hsBa where 
Ba — 2D /[is is the maximum uniaxial anisotropy field. E.g., 
for FesoPtso we have TP rcc = 5 ps, the maximum anisotropy 
field is ~ 7 T and the magnetic moment per nanoparticle is 
around 22000/iB [26]. The field-free solution of (fTJ is known; 
e.g. for a uniaxial anisotropy and starting from the angles 
(j)f{t = to) and 9 fit = to) one finds (e.g., 13011 ) 
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t-to 
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1 + tan 2 9 f (t )-e 



2a(t-t ) 



1 + COS 

7X7^ (*-*°) 



tan#f(t) = tan 0/ (to) • e 

(2) 

"+" ("-") refers to < 9 < tt/2 ( tt/2 < 9 < n). 
To control the dynamics we apply along the x and y axis two 
magnetic field pulses b x and b y of durations 2e and shapes 
/(f) centered at some moment f = fo- Their relative strengths 
is given by the mock angle <po, with tan^o = M/IM^ the 
total fields strength is |/|6n/(2e). Hence bo(t) = h x + h y is 



2 



; ' ' b ° (cos o e x + sin 0oe y ), t - e < t < t + e 



bo(*) H 2e 



Switching to a new time variable r(t) 
rive for the equation of motion 



elsewhere. 



t-(t +e)+2e 
2e 



(3) 

we de- 
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+ a cos cos i 



where 6<fi = <fr — <po and p = 1/(1 + a 2 ). If the magnetic 
pulses are shorter than the precessional period then from eq. 
© we infer for the angles stroboscopic evolution from before 
(</>(i~), 6(t~)) to after (4>{t + ), 9(t + )) the pulses the relation 
(we introduced t~ := to — £, t + := to + e) 
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a cos cos ( 



which is valid up to terms of the order (e/T prec ) 2 . After the 
pulse, i.e. for t > t + the dynamics is governed by eq.© 
with the initial conditions (j>t = <f>(t + ), Of = 9(t + ). This 
procedure is repeated by applying further pulses accordingly. 

Controlled switching.- As we are interested in switching 
we require in the spirit of local control theory that 



0(t+)>9(t-) vt+r 



(6) 



procedure is repeated until we achieve the state with = ir/2. 
As clear from © the tilt angle is always increased upon the 
pulse with an amount that goes linearly with the fields strength 
bo. On the other hand the variation of <f> with bo is only loga- 
rithmic, in fact if the time delay between the pulses is only a 
fraction of the precessional period, tj> is hardly changed. 

Freezing. - The present LCT scheme allows also for the sta- 
bilization of the magnetization around a desirable target non- 
equilibrium angle 9 t : At first, starting from a given state 
we apply the control scheme and achieve 9 t at some time 
t t . During a field-free period r the angle 9 t develops to 
f(tt + t). To compensate for this change we apply a pulse 
(centered at £o,t) according to our scheme this shifts the angle 

to 6+ = f (t t + t) + fcn /+a^ ) - To stabilize the magnetiza- 
tion we choose bo such that 9 + = 9 t . The procedure is then 
repeated during the stabilization time. To minimize the ad- 
justment of bo between consequent pulses the repetition rate 
should be large. 




FIG. 1: (Color online) Evolution of 9(t) according to the proposed 
control scheme and for <j>(t = 0) = 7r/180 = 9(t — 0),(j>o = 
a,Tctw(b y /b x ) = 27r/3, a = 0.05, / = 1, bo = 0.2. Inset shows 
the short-time behaviour (pulses are off for 9 > 7r/2). 



As inferred from eq.©, this condition is fulfilled if 5(f) = <f> — 
4>o = 37r/2. If a sequence of the pulses © each centered at 
the times to,i is applied then S(t) evolves as 



tan i — h 2 1+Q 2 



tan 



baf{t 0j i) 
1 + a 2 ' 



(7) 



where tf — ta,% ± £■ 

The realization of this LCT scheme is then as follows: Starting 
from a known (e.g. equilibrium) state <f> — 0(0); — 9(0) we 
apply at t = <o,i me m "st fields b x and b y © with strengths 
such that (f>o = (f)(0) - 3vr/2 (cf. Fig.©). Eq. © deliv- 
ers the tilt angles 6(tf) and <j>(t^). During a time lag (dark 
time) T\ the propagation proceeds according to eq.© with 
the initial values 4>f(i ) = <j>(tf) and 0/(*o) = Q(tf). At 
t = to, 2 we apply a second pulse with b x and b y such that 
(f>o = <f>f{ti + t\) — 3ir/2. From eq. © we deduce that af- 
ter the second pulse 9(t+) = 9 f (t+ + n) + b °/+^ 2) ■ This 
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FIG. 2: (Color online) #(t) for different pulse durations (solid rect- 
angles). T prcc is the precessional period and bo = 0.3, a = 0.05. 



Numerical results and illustrations.- Fig. [T]shows the magne- 
tization reversal according to our zero temperatures (T = 0) 



3 



analytical scheme and in the damping regime appropriate for 
magnetic nanoparticles. Fig. [T]confirms our analysis and the 
physical picture drawn above. However, the following issues 
need to be clarified for this procedure to be of practical in- 
terest. 1.) Do we need a precise tuning of the pulses du- 
rations, 2.) will thermal fluctuations invalidate our findings, 
and 3.) how effective is this scheme when applied to other 
type of anisotropy fields. To address these points we imple- 
mented a finite-temperature full numerical realization 13111 of 
the present control scheme (cf. |H Q] and references therein 
for an overview on numerical micromagnetic methods), i.e. 
the analytical expressions deliver the appropriate input param- 
eters for the numerics. The damping parameter is chosen ac- 
cording to experimental findings [2]. For the simulation pre- 
sented here we use square-shaped pulses, i.e. f(t) = 1 for 
to — e < t < to + e. Basically the same conclusions are 
valid for other pulse shapes, e.g. Gaussian pulses ll32Tl . Fig. 
|2] demonstrates the evolution sensitivity of the angle 9 when 
pulses with different durations are applied. It also shows the 
range of validity of our scheme. As inferred from Fig. |2] 
a fine tuning of the pulse duration is not mandatory as long 
as it is smaller that T prcc . The strength bo determines the 
value of the tilt angle (as follows from eq.©). The insensi- 
tivity to the pulse duration is favorable for practical applica- 
tions, however the generation of magnetic pulses shorter than 
T prcc might be a challenge; the light-induced generation of 
sub-picosecond shaped magnetic pulses [33] may circumvent 
this problem. As for the role of the magnetization dynamics 
during the pulses our full numerical simulations (cf. Fig. [3) 
confirm qualitatively the analytical predictions. According to 
e q-dZl* a minimal fields strength bo is required for switching, 
for bo determines 0(t + ). To realize the stabilization scheme 
outlined above one tunes bo to steer the magnetization to a 
non-equilibrium 9 t (cf. Fig. [3) and keep it there (as long as 
b is on). 
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FIG. 3: (Color online) Tilt angle 8(t) within the present the local 
control scheme for different fields strengths bo. Other parameters: 
a = 0.05, T = OK, (Pulses are off when 6 > tt/2). 



pulses have to be applied even if t > tt/2, since due to ther- 
mal excitations the magnetization may swing back to the orig- 
inal state. This effect is avoided by applying the pulses even 
if > tt/2 (Fig. 2J lower panel). Generally, we observe that 
thermal fluctuations have little influence on the effect of the 
pulses (i.e., on the dynamics during and right after the pulses), 
in contrast to continuous fields 13111 . On the other hand the 
field-free processional motion between the pulses is generally 
modified at T > 0. 

The possibility of field-assisted stabilization (freezing) can be 
exploited for fast field-assisted thermal switching: Starting at 
T w we utilize our scheme to drive the magnetization to a 
state t < 7r/2 (as shown in Fig. |5]l and then freeze it there. 
At low temperatures switching does not occur irrespective of 
the waiting time (inset of Fig. [5]). When the temperature in- 
creases however, the thermal fluctuations increase but can not 
lead to a reversal in absence of the field, as demonstrated by 
the inset of Fig. [5] The presence of the fields assists a fast 
magnetization reversal, a behaviour that can not be realized 
with static fields, since a magnetization freezing is necessary. 
In practice, the reversal process may be functionalized a fast 
thermal sensor to monitor swiftly a temperature increase. 
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FIG. 4: (Color online) Temperature-dependent controlled evolution 
of the angle 6(t) (a — 0.05, bo — 14.77). The pulses are applied if 
9 < tt/2 only (top panel) or throughout (below). 
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Figure@]proves the robustness of the scheme to thermal fluctu- 
ations. Here we highlight a special feature of the temperature- 
dependent magnetization dynamics: To achieve switching, the 



FIG. 5: (Color online) Thermal-assisted controlled switching in the 
presence of short pulses with an amplitude bo = 8.86. Inset shows 
switching is not possible for bo = 0. 
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The question of to what extent the present scheme is appli- 
cable to another anisotropy type we address by studying the 
magnetization control of Fe7QPt30-nanoparticles which pos- 
sesses cubic anisotropy 11271 13411 . For a cubic anisotropy the 
field-free ground-state energy landscape contains several min- 
ima 113511 . By switching we mean then a magnetization trans- 
fer between these minima and not necessarily a change from a 
parallel to an antiparallel state, as in the uniaxial case. Fig. [6] 
demonstrates the applicability of our control proposal. Start- 
ing from a state close to one energy minimum the magne- 
tization precesses and relaxes in a field-free manner to the 
ground state. When the magnetic pulse is applied according 
to our LCT the magnetization is transferred almost directly to 
the next energy minimum in the positive energy semi-sphere. 
With the freezing scheme outlined above it is even possible 
to stabilize the magnetization on top of the barrier (Fig. |6). 
Summary.- A sequence of two perpendicular magnetic pulses, 
each with a duration less than the precessional period is ca- 
pable of increasing monotonically the magnetization tilt angle 
as to achieve a predefined state. As shown analytically (for 
T = 0), this is possible if the relative strengths is tuned ap- 
propriately. Full numerical simulations accounting for finite 
temperatures and different types of anisotropy fields demon- 
strated the usefulness of this scheme for magnetization "freez- 
ing" and switching within tens of picoseconds. We also illus- 
trated how the method can be exploited for picosecond field- 
assisted thermal-switching. 




FIG. 6: (Color online) Polar diagram of the energy surface for a cubic 
anisotropy with magnetization trajectories. Left panel is a top view 
on the energy surface: For bo = (dark trajectory); For a bo = 2.06 
control field (light trajectory). Trajectories start at cj>(t = 0) = 1.97T, 
6(t = 0) = 7r/3.8. Right panel is a bottom view at the energy 
surface: Freezing field is bo = 0.59 and the magnetization is initially 
at the position marked (X). In both cases a = 0.05. 
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